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Abstract. Wetting phenomena, molecular protrusions of lipid bilayers and membrane stacks under 
lateral tension provide physical examples for interacting surfaces with tension. Such surfaces are 
studied theoretically using functional renormalization and Monte Carlo simulations. The critical 
behavior arising from thermally-excited shape fluctuations is determined both for global quantities 
such as the mean separation of these surfaces and for local quantities such as the probabilities for 
local contacts. 
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1 Introduction 

The surfaces considered here are interfaces and membranes 
which are governed by tension. An interface or domain 
wall which represents the stable contact region between 
two bulk phases of matter is always characterized by a fi- 
nite interfacial tension, see, e.g., Q. Flexible membranes, 
on the other hand, are sheets of amphophilic molecules 
such as lipid bilayers and are usually controlled by cur- 
vature and bending rigidity, see, e.g., 0. However, these 
membranes may be subject to an applied lateral tension 
arising from external forces or constraints. In addition, 
lipid bilayers exhibit molecular protrusions on small scales 
which are also governed by an effective tension. 

In this paper, we study the interactions of two or sev- 
eral such surfaces and focus on the renormalization of 
these interactions arising from thermally-excited shape 
fluctuations. Physical examples are provided by (i) Wet- 
ting layers which are bounded by two interfaces; for re- 
views, see In this case, the interactions between 
these two interfaces are renormalized by capillary waves 
j?J ; (ii) Lipid bilayers at separations which are small com- 
pared to the bilayer thickness. Such a situation arises for 
stacks of bilayers in the presence of a relatively large ex- 
ternal pressure acting on these stacks. In this case, the in- 
teractions between adjacent bilayers are renormalized by 
molecular protrusions - thermally excited displacements 
of neighbouring lipid molecules which change the area of 
the lipid- water interface ||; and (iii) Swollen bunches or 
stacks of membranes at separations which are large com- 
pared to the thickness of these membranes. This corre- 
sponds to the limit of relatively small external pressures 
acting on the membrane bunch or stack. Previous theo- 



retical work on such bunches has focused on tensionless 
membranes as reviewed in Ref. ||. In contrast, we will be 
concerned here with membranes under lateral tension. 

The paper is organized as follows. In Section |^, we will 
define the theoretical models appropriate for (i) complete 
wetting, (ii) protrusion forces between bilayers, and (iii) 
membranes under lateral tension and explain that these 
models are equivalent to each other if the various param- 
eters are identified in an appropriate way. In Section |^, 
we define the global and local critical effects which we 
will study in the remainder of the paper. One critical ef- 
fect which has not been studied previously in a systematic 
way is the behavior of the probability for local contacts be- 
tween the surfaces. Our results for two surfaces with hard 
wall and soft wall interactions are described in Sections ^ 
and ||, respectively. 

In order to simplify the presentation, we will explicitly 
discuss our results in Sections [| and || using the terminol- 
ogy of molecular protrusions for lipid bilayers. However, 
if one makes the appropriate identification of parameters 
as explained in Section @, these results can be directly ap- 
plied both to complete wetting and to membranes under 
lateral tension. 



2 Theoretical framework 
2.1 Wetting phenomena 

Consider an interface between two bulk phases, say a and 
7, and let us change a thermodynamic field, such as e.g. 
temperature or pressure, in order to move the system to- 
wards a triple point where a third phase, /?, can coexist 
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with the two phases a and 7. Then, in thermal equilib- 
rium, a layer of the (3 phase may appear in the a(3 inter- 
face. As one comes closer and closer to the triple point, 
the thickness of this intermediate layer may continuously 
grow: this is the case of complete wetting. On the other 
hand, no layer may appear in the a(3 interface or the 
thickness of this layer may saturate and remain finite as 
the triple point is attained: this is the case of incomplete 
wetting. 

When the (cry) interface contains a wetting layer of f3 
phase, it splits up into an (a/3) and a (Pj) interface. The 
mean separation, I > 0, of these two interfaces is equal to 
the thickness of the wetting layer. Likewise, the excess free 
energy of the layer can be regarded as an effective inter- 
action between these interfaces. Here, we will focus on the 
situation where the interaction potential between these 
two interfaces is purely repulsive and short-ranged. An 
example is provided by complete wetting in 3-dimensional 
lattice gas or Ising models. 

2.1.1 Solid-on-solid models for wetting 

Lattice gas or Ising models on a cubic lattice with a pla- 
nar surface presumably provide the simplest models for 
wetting in three dimensions even though the correspond- 
ing phase diagram is already quite complex, see, e.g., pcj| . 
Wetting occurs at or close to two-phase coexistence pro- 
vided the bulk of the system is in one phase, say in the 
spin-up phase (corresponding to the dilute vapor phase 
in the lattice gas) whereas the surface of the semi-infinite 
system favors the other phase (corresponding to the dense 
liquid phase in the lattice gas). In such a situation, one has 
a thin wetting layer of spin-down phase which is bounded 
(i) by the interface between this layer and the bulk phase 
and (ii) by the planar surface of the semi-infinite system. 

At low temperatures, the interface is essentially flat. 
As one increases the temperature, typical excitations con- 
sist of islands on this flat interface. These excitations can 
be studied in the framework of so-called solid-on-solid 
models in which the planar surface of the semi-infinite 
system is described by a square lattice with lattice sites 
Xi and lattice constant a, and the local separation of the 
interface from this surface is described by the separation 
field li = l(xi). The configurational energy or effective 
Hamiltonian of this separation field is given by 

H{1} = $>J/a)|Z, - lj\ + XX (2.1) 

m 

where (ij) indicates a summation over all nearest neigh- 
bors, J is the nearest neighbor coupling constant in the 
Ising model, and V(l) denotes the effective potential act- 
ing on the interface. |ll| For complete wetting, this po- 
tential has the simple form 

V(l) = V hw (l)+Hl (2.2) 

where the hard-wall potential V/ lw (0 is defined by V/ lw (0 = 
for I > and Vh w (l) = 00 for Z < and H is the bulk 



magnetic field in the original Ising model. At temperature 
T (measured in energy units), the statistical weight for the 
configuration Z, = l(xi) is given by the Boltzmann factor 
~ exp[-W{Z}/T]. 

2.1.2 GauBian interface models for wetting 

The model as given by (|]l]) an d ( ] 2.2| ) has been studied 
by Monte Carlo simulations. fl^ , |l3| In these simulations, 
the mean separation of the two surfaces was determined in 
the limit of small magnetic field H. It was found that this 
quantity exhibits the same critical behavior in the discrete 
solid-on-solid model as in the corresponding GauBian in- 
terface model. In its discrete version, the latter model is 
defined by the effective Hamiltonian 

mi} = J2 \m h? + X « 2 ^) ( 2 - 3 ) 

(ij) ' 

where S represents the interfacial tension (or stiffness) of 
the fluctuating interface on large scales. In the continuum 
limit, this leads to 

H{1} = J d 2 J^(VZ) 2 + F(Z)} (2.4) 

where both x and I are continuously varying coordinates. 
The latter model has been studied by functional renor- 
malization as discussed in Ref. ||. 

The critical behavior of the mean separation (l) is iden- 
tical in the solid-on-solid and in the Gaufiian model pro- 
vided one identifies the interfacial tension with 

S = c s (2J/af/T with c s ~ 2.6 (2.5) 

for the square lattice. 

In the present work, we will address the question to 
what extent this correspondence between the two models 
is also valid for other quantities such as, e.g., the proba- 
bility for local contacts. In fact, we will see that the latter 
quantity has the same behavior for hard wall interactions 
but exhibits different behavior for soft wall interactions, 
see Sect. 5 below. 



2.2 Molecular protrusions of bilayers 

Now, let us consider a rather different system consisting 
of several lipid bilayers. It has been well-established by 
many experiments that such bilayers in aqueous solution 
experience strong repulsive forces at small separations of 
the order of 1 nm. jf4| It was originally thought that this 
short-ranged repulsion represents a hydration effect and 
reflects the perturbed water structure in front of the polar 
head group More recently, it has been shown, how- 
ever, that a similar repulsive force can also arise from the 
molecular roughness of the lipid-water interfaces. 

Lipid bilayers consist of rod-like molecules with a hy- 
drophilic head group and usually two hydrophobic tails. 



Andreas Volmer et al.: Critical behavior of interacting surfaces with tension 



3 



These rods are packed in such a way that the hydrophilic 
head groups shield the tails from the surrounding wa- 
ter. The corresponding lipid/water interface is roughened 
by thermally-excited fluctuations as has been observed in 
computer simulations |l6|,|l7| and has been deduced from 
scattering experiments J18| , |19| , p0[ . These thermal fluctua- 
tions correspond to protrusions, i.e, to relative displace- 
ments or deformations of the lipid head groups. Since these 
protrusions will, in general, change the surface area of the 
lipid-water interface, they are governed by an effective 
tension. 

The repulsive force arising from the protrusions of sin- 
gle molecules was calculated in Ref. pi] . The effect of 
collective protrusions which involve whole groups of lipid 
molecules was first studied in Ref. Q within the frame- 
work of discrete and continuum interface models. 



2.2.1 Models for collective protrusions 

Consider one such lipid/water interface formed by the 
head groups of the lipid molecules which is located in 
front of a planar surface or wall. The head group at po- 
sition Xi has the local separation li from the wall. If the 
lipid/water interface is flat, its separation from the wall 
is constant and all U are equal. Thermally-excited pro- 
trusions roughen this flat interface and lead to position- 
dependent configurations U = l(xi). 

The configurational energy of these protrusions con- 
sists of two contributions. The first contribution is given 
by the excess free energy arising from the relative displace- 
ments of the molecules which is governed by the interfacial 
tension S . The latter quantity represents the free energy 
per unit area for the interface between the nonpolar part 
of the molecule and the water. In the discrete model to 
be defined, each molecule is supposed to have n nn near- 
est neighbors. The molecular shape is approximated by 
a rigid column with circumference a and cross-sectional 
area A a . For n nn = 6, one has hexagonal columns on a 
triangular lattice with lattice constant a = a /2\/3 and 
cross-sectional area A Q = V3a 2 /2. 

The second contribution to the configurational energy 
of the protrusions represents the excess free energy aris- 
ing from the interaction of the lipid/water interface with 
the planar surface or wall. The corresponding interaction 
potential V(l) has the generic form 



V(l) = V hw {l) + V hy exp[-l/l hy ] + PI 



(2.6) 



The first term of V(l) is again the hard wall repulsion 
which ensures that the two s urfa ces cannot penetrate each 
other. The second term in ( |2.6| ) represents the hydration 
interaction arising from the perturbed water structure in 
front of the lipid head groups; it is parameterized by a 
certain decay length ly ty and a certain amplitude Vhy The 
last term PI is the excess free energy arising from the 
external pressure P. For finite P, the lipid/water interface 
will have a finite mean distance from the wall and will 
undergo an unbinding transition in the limit of small P. 



If one combines both contributions, one obtains the 
configurational energy or effective Hamiltonian |s) 

H{1} = J2 (a S /n nn ) \l t A oV(k) (2.7) 

<ij> i 

for the collective protrusions described by the configura- 
tion li = l{xi). Comparison with the solid-on-solid model 



(2.1) shows that both models are completely equivalent. 

The model as given by (2/7) was studied by Monte 
Carlo simulations for a triangular lattice with n nn = 6; in 
addition, the corresponding Gaufiian model was studied 
by functional renormalization. Q It was again found that 
the critical behavior of the mean separation is the same in 
both models provided one chooses the effective interfacial 
tension 



Aj = En 



c s (a S /n nn ) jT. 



(2.8) 



For the triangular lattice with n nn = 6 nearest neighbors, 
the Monte Carlo simulations in Refs. || lead to the es- 
timate eg ~ 2.4 which is different from but close to the 
value c-s — 2.6 as obtained for the square lattice. 



2.3 Membranes under lateral tension 

On sufficiently large scales, flexible membranes which do 
not experience any external force or constraint are gov- 
erned by their bending rigidities and curvature energies, 
see the various reviews in Ref. 0. On the other hand, 
the presence of such forces or constraints tends to induce 
lateral tensions which reduces the bending undulations 
of the membranes ]22] , ^3| . One example is provided by 
a closed vesicle which is inflated by an osmotic pressure 
difference AP = Pi n — P ex . The surface of this vesicle will 
attain an essentially spherical shape with radius R ve and 
the corresponding tension S satisfies the Laplace equation 
S = R ve AP/2. 

If the vesicle is unilamellar, i.e., if it is bounded by a 
single membrane, this membrane is subject to the tension 
S. If the vesicle is multi-lamellar, i.e., if the vesicle surface 
consists of N closely packed membranes, each membrane 
should experience a lateral tension of the order of S /N . 

It is also possible to apply a controlled suction pressure 
AP to the vesicle by micropipet aspiration 0,^5) . In this 
way, one can directly control the lateral tension S. For a 
multi-lamellar vesicle, one then has N interacting mem- 
branes under lateral tension. In addition, the tense vesicle 
can also be pushed, via the micropipet, against a planar 
substrate or wall. Within the contact region, one then has 
an oriented stack interacting with the planar surface. 



2.3.1 Continuum GauBian models 

Now, consider an oriented stack of membranes labeled by 
n with 1 < n < N . Each membrane is characterized by 
its bending rigidity n n and is taken to experience the lat- 
eral tension E n . For an oriented stack of membranes, the 
position of membrane n can be described by the height 
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variables h n {x) which measures the distance from a ref- 
erence plane with coordinate x. The interaction potential 
between two neighboring membranes is again denoted by 
V(hn+i — h n ). Using this parameterization, the effective 
Hamiltonian for the whole stack has the form 



H{h} = j d 2 x [ J2 

J ^ n=l 



2 

JV-l 



1 s n (v/i„) 2 + (y 2 h 7 , 



XV — J. x 

V(h n+1 - h n ) I 

n=l > 



(2.9) 



For N — 2, one has a rigidity-dominated regime for suf- 
ficiently small scales and a tension-dominated regime for 
sufficiently large scales [pj2(|. In the latter regime, the 

bending terms ~ (V 2 ^„) become irrelevant. Such a be- 
havior is to be expected for general N > 2, and the crit- 
ical behavior is then governed by the interplay between 
the tension terms and the interaction terms. In order to 
study this interplay, it is convenient to make an orthogo- 
nal transformation from rescaled height variables y/ S n h n 
to new fields z n which contain the 'center-of-mass' coor- 
dinate 



ZN 



N 

E 



N 



n=l 



1/2 



(2.10) 



The latter coordinate decouples from the other fields z n 
with n < N since it does not enter in the interaction terms 
~ V(h n+ i — h n ). In this way, one arrives at models for the 
N—1 fields z n with n < N which are linear combinations 
of the N—1 separation fields l n = h n +i — h n . For N = 2, 
one obtains a model for the single field z\ ~ h-2 — hi which 
is equivalent to t he continuum Gaufiian model for wetting 
as given by (2.4). 



3 Critical behavior 

In the following, we will discuss our results using the ter- 
minology which is appropriate for molecular protrusions 
of bilayers. The potentials V(l) considered here consist 
of repulsive hard-wall and soft -wall interactio ns b alanced 
by an external pressure term PI as given in (2.6). In the 
present section, we will define the global and local quanti- 
ties which become critical as the pressure P goes to zero. 



3.1 Global quantities 

Thus, consider two surfaces governed by tension which are 
pushed together by the external pressure P. The config- 
uration of these two surfaces is described by their local 
separation I = l(x). As P is decreased to zero, the two 
surfaces unbind which leads to the divergence of several 
length scales. First of all, the parallel correlation length 
£|l diverges as 

£|| ~ l l p ^ (3-1) 
which defines the critical exponent uu . 



In addition, both the mean separation I = (I) of the 
two surfaces and the roughness 



^ = <(Z-(i)) 2 > 



1/2 



(3.2) 



of the separation field I are singular for small P. For 
protrusions, w hich are governed by the effective tension 
£ = S pr as in (|2.8| ) , the interfacial roughness is related to 
f|| via 

(T/2nS)^yjH^/a) (3.3) 
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where a denotes the small-scale cutoff as before . This 
relation is analogous to the well-known scaling of the 
roughness of a free interface as a function of its linear 
dimension. Note that the basic length scale which sets the 
size of the interfacial roughness is given by the protrusion 
length 



It. 



(T/2ttU) 1/2 



(3.4) 



For three surfaces, one has two separation fields l± = 
h-i — hi and I2 — ha — h%. In this case, one convenient 
choice is to use the two fields di ~ h + h, which is the 
thickness of the whole bunch, and d,2 ~ li — li- 



3.2 Local contacts 

As the surfaces unbind, the probabilities for local contacts 
of the surfaces decay to zero. More precisely, the proba- 
bility 7^26 for pair contacts is found to behave as 



V 2b - 1/& - P" 2 



(3.5) 



which defines the critical exponents v-i- This exponent sat- 
isfies the scaling relation 



V 2 = C,2V 



(3.6) 



For a bunch of three or more interacting surfaces, one 
may consider the probabilities V n b for local n-contacts. 
These quantities have the scaling behavior 



Vnb 



with 



0|| 



(3.7) 



Local contacts have also been studied for 1-dimensional 
lines (or strings or directed walks) governed by line ten- 
sion |B7 28,29LpO| and for tensionless membranes [B0|. In 



the following sections, we will study these quantities for 
two and three interacting surfaces governed by tension. 



4 Two surfaces with hard wall interactions 

In this section, we study the critical behavior of two sur- 
faces interacting with hard wall interactions Vh w (l) with 
Vhw(l) = for I > and VhwQ) = 00 for I < as before. 
Since this interaction contains no length and no energy 
scale, it does not introduce any dimensionfull parameter 
into the model. 
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We begin with a brief summary of the results from a 
functional renormalization group treatment, which gives 
predictions for the mean distance I and the roughness £j_ 
of the interfaces. Additionally, we consider the pair con- 
tact probability V 2 b and derive predictions for its scaling 
behavior. Monte Carlo (MC) simulations are used to test 
the theoretical results. The GauBian and the solid-on- 
solid model are treated separately. 



4.1 GauBian models 



First, we give a short review of the functional renormaliza- 
tion group treatment of the GauBian model as defined by 
(2.4), following The functional renormalization group 
is used to integrate out thermal fluctuations on scales be- 
tween the microscopic length a and the lateral correla- 
tion length £|| , which will be finite for finite values of the 
pressure P. Technically, this integration is performed in a 
one-step nonlinear functional renormalization group with 
a rescaling factor b = £||/a, using a formalism which is an 
extension of Wilson's approximate recursion relation |]3lf , 
and which has been developed in studies of the wetting 
transition Q. This procedure yields an effective interac- 
tion potential on the scale of the lateral and perpendicular 
correlation lengths £|| and Now, on this scale, mean 
field theory can be used, allowing to set up a self consis- 
tent procedure to determine the mean separation I = (/) , 
and £n by minimizing the effective potential. This effec- 
tive potential is now, on the scale £|| , defined as a simple 
superposition of the effective hard wall potential and a 
linear pressure term: 



v CB {i)^vZ{i) + Pi, 

where V^(l) = — v lnerf(Z/\/2^j_). The parameter 
T/A a represents the basic energy scale. 
Minimization with respect to I leads to 



(4.1) 



P 



01 



i=e 



V2v 

7^1 



-r/2?i 



(4.2) 



In addition, we have the mean field relation between £|| 
and the curvature of the effective potential in the mini- 
mum position 

S 



d 2 v, cS 

hw 



or 2 



(4.3) 



and the rel ation between the two length scales £n and £j_ 
as given in (3.3). Solution of this set of equations yields a 
relation between the pressure and the mean separation I 
that is given by 



P ~ Phw& ^ P (Jpr I ^ 7 



(4.4) 



with the protrusion length l pr — (T/2t:E) 1 ^ 2 as before. 
Likewise, the surface roughness £j_ is found to obey the 
relation 



P ~ Phw±e 



The two pressure amplitudes take the form 



Phw = 



rV4, 



and 



hwA- 



(4.6) 



Inverting the two relations (4.4) and (4.5), one finds 

I « l pr [ln(P hw /P) - 1/4 lnln(P,™/P) ] (4.7) 

and 



ei « i 2 pr /2 HP hl 



jp) 



(4.8) 



in the limit of small P. 

The mean separation and the correlations lengths rep- 
resent global quantities of the interface. In order to intro- 
duce a local quantity, namely the pair contact probability 
P2b, we define the probability distribution 



T A (l) = (9(l + Al)-8(l))/Al, 



(4.9) 



where 0(1) is the Heaviside step function and Al is a mi- 
croscopic length scale. Va{1) gives the probability to find 
the surface height l(x), at a given position x, in the inter- 
val I . . . I + Al. Then, the probability V^b for pair contacts 
is given by 



V2b=VA(l = 0). 



(4.10) 



In the following, we will omit the subscript A, and take Al 
to be a fixed microscopic length scale Al. If one assumes 
that V(l) is well approximated by a GauBian distribution 
which is centered at the mean separation I and has the 
width ~ the pair contact probability is given by 



V 2b - ex P M 2 /2£i]/6 



(4.11) 



Note that in the GauBian case, the full distribution "P(Z) 
is fully determined by its first two moments, £ and y.. 
Inserting the asymptotic p ress ure de pend ence o f thes e two 



into (1.11), one 



length scales, as given by (4.7) and (4 
immediately obtains 



V 2b ~ exp [- ln 2 (P^/P) / ln(P hwX /P)] (4.12) 
- P V2 with v 2 = 1, (4.13) 

neglecting confluent logarithmic corrections to the leading 
scaling behavior. The contact exponent v 2 , which has been 
introduced in Sect. 3.2, governs the asymptotic behavior 
for P > P hw and P > P hw± . 

The above assumption about the GauBian character 
of the probability distribution V(l) has been confirmed by 
Monte Carlo simulations of the discrete GauBian model. 
In the simulations, the discrete form 



H{z} 



H{z} ^ ( Zi - Zj f 



E 



2n r , 



/XUhwjZj) + QZi) 



(4.14) 

of the effective Hamiltonian has been studied, where we 
have used the dimensionless quantities z = l/l sc and Q = 
P/P sc with 



(4.5) l sc = ^T/Sn 

r> 



and 



P,, 



^TEn nn /a 2 . (4.15) 
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Fig. 4.1. Some typical values for the relative error of the 
measured mean distance as a function of the simulation time, 
measured in MC steps per site At. The SOS model with one 
flexible surface fluctuating against a hard wall was simulated, 
with a system of linear dimension Lm = 64. 



Fig. 4.2. Probability distribution V(z) for the discrete 
Gaufiian model of a fluctuating surface which is pushed against 
a hard wall by the external pressure Q. The data points are 
fitted by Gaufiians. 



For all MC simulations, a triangular lattice with 6 
nearest neighbours per site was used, and the simulated 
area had the form of a parallelogram. Periodic bound- 
ary conditions were applied in the direction parallel to 
the base line of the parallelogram and in the direction 
perpendicular to this line. The linear dimensions in both 
directions were always taken to be equal. Typical lattice 
sizes were Lm=32 and 64, where L\\ is the linear dimension. 
For optimization of calculation time, the height variables 
Zi were discretized; we have checked, however, that the 
results do not differ from those obtained when using float- 
ing point variables. Finally, the acceptance rate of the MC 
step was continuously controlled by tuning the maximum 
height variation per step and site. The length of the MC 
runs varied between 10 5 and 10 6 MC steps per site, de- 
pending on the chosen value of the external pressure Q. 

The statistical errors of the MC results were controlled 
in the following way. Starting with a small MC sampling 
interval At, the mean values of all observables were mea- 
sured within n equal intervals of length At/n, giving n 
values with /j, = 1, . . . , n. Usually, we have chosen the 
value n = 10. From these values, the statistical error 5A, 
defined as 8 A 2 = ( [ ± £™ =1 (A^~(A))} 2 ), was calculated. 
Then, At was multiplied by n, and the procedure was re- 
peated on the larger time scale. In this way, one can plot 
the relative statistical error S Aj (A) of the observable A as 
a function of At. The results of this procedure are illus- 
trated in Fig. 4.1. Note that when At becomes larger than 
the relaxation time, the relative error decreases monoton- 
ically. In the MC results presented in this paper, we have 
plotted no error bars because the statistical errors were 
typically smaller than the data symbols. 

MC data of probability distributions 1^(1) fo r th ree dif- 
ferent values of the pressure Q are shown in Fig. 4.2. These 
data were obtained in simulations of one flexible surface 
pushed against a hard wall, and are well fitted by Gaufiian 
distributions. 

Consequently, the value v 2 = 1 of the contact exponent 
is confirmed by MC simulations of the discrete Gaufiian 



model. Similar data have been obtained in simulations of 
the solid-on-solid model as explained in the next section. 

In addition, the p redi ctions for the global quantities £ 
and as given by (47) and (T_8), have been confirmed 
in our simulations. The corresponding data are not shown 
here since they are similar to the corresponding data for 
the solid-on-solid model shown in the next section. From 
the fits of the mean distance £ to the functional form as 



given by (4-7), we have extracted the dimensionless pro- 
trusion length z pr = lpr/l sc = 0.55±0.02. We have checked 
that this value of consistent with the fits of the 



roughness £j_ as in (4.8) 



4.2 Solid-on-solid models 

In order to set up Monte Carlo simulations of the solid-on- 
solid model for collective protrusions as defined by (2.7), 
we use the dimensionless length z = l/l sc and pressure 
Q = P/P sc , where the length and the pressure scale are 
now given by 



l sc = T/a E and P sc = a E /A . (4.16) 

In terms of these rescaled quantities, the effective Hamil- 
tonian for one protruding surface reads 

H{z} =^2\ Zi - Zj\/n nn + ^(Qzi + U hw (zi)), (4.17) 

(ij) * 

where n nn is the number of nearest neighbours of each 
molecule as before. The generalization for two protruding 
surfaces is obvious and leads to 



T-i{zi, z 2 } = — !— (|zi,; - zx.jl + 

^ — ' n 



'•2. 



«2,il) 



■ >^ Q - z 2 ,i) + U hw {z lti - z 2 ,i), (4.18) 
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Fig. 4.3. MC data for the SOS model. Mean distance I (left) 
and roughness £x (right) versus pressure Q, both in the case 
of one and of two fluctuating surfaces. The mean distance is 
fitted by £ ~ V(ln(Q h „/Q)-l/41nln(Q h „/Q)), see Eq. @, 
with i pr = 1.57 and Qhw = 0.47 for one surface and l pr ~ 1.06, 
Qhm = 0.81 for two surfaces. The values fo r th e roughness are 
fitted by ~ l 2 pr /2 ln(Q hw± /Q), see Eq. Q, with Q hw± = 



0.74 for one surface and QhwX = 1-55 for two surfaces. 
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Fig. 4.4. Mean distance (z) and roughness £x versus pressure 
Q in the high pressure regime for the SOS-model. The MC data 
were obtained for a pure hard wall potential and a triangular 
lattice with linear dimension L« = 32. 



which, after introducing the difference coordinate z = Z\ — 
zi and the center-of-mass coordinate z = z\ + zi 1 takes 
the form 

H{z, z} — ^2 —— Max(|zi - zj\, \xa - Zj\) 



Zi + U hw (zi). 



(4.19) 



Note that in the solid-on-solid (SOS) model, the center- 
of-mass coordinate and the difference coordinate do not 
decouple, in contrast to the Gaufiian model. By this decou- 
pling, the Gaufiian model with two surfaces can be exactly 
mapped to a model for two independent surfaces, one of 
which is free (the center-of-mass coordinate), the other 
one being subject to the external potential (the difference 
coordinate). We have employed Monte Carlo simulations 
in order to check whether the critical behavior of global 
qua ntitie s in the SOS model for two interfaces as given 
by ( |4. 18| ) , is nevertheless consistent with the behavior ob- 
tained from the Gaufiian model. 

This assumption is indeed confirmed by the results 
from MC simulations, as shown in Fig. [O. There, the 
mean distance £ and the roughness £j_ are plotted versus 
the pressure Q, both in the case of one and of two fluc- 
tuating surfaces. The length scale l pr (2) in the latt er c ase 
relates to the protrusion length l pr as defined in ( |3.4[ ) in 
the following way, using the terminology of the membrane 
model: Consider the surface tension S of a single mem- 
brane surface in the Gaufiian model. First, in the case of 
two identical surfaces, the reduced surface tension S T of 
the distance coordinate l\ — I2 is given by S/2. Second, 
the microscopic tension Sq has to be multiplied by 2 since 
two adjacent molecules in the lipid bilayer should move 
coherently in order to prevent the generation of energeti- 
cally unfavorable cavities between them. With £ <~ Sq as 
in (fO|), these two effects lead to l pr (2) = l pr /V^- 

The actual values for z pr , as extracted from the simula- 
tions results in the same way as described in the preceding 
chapter on the Gaufiian model, = 1.57 ±0.03 for 



one surface and z pr ( 2 ) — 1.07 ± 0.02 for two surfaces, 
corresponding fits are shown in Fig. 



The 



4.3 



Let us shortly focus on the high pressure regime with 
Q > 1. In this regime, the dimensionless external pressure 
Q is larger than the microscopic surface tension (which 
has been scaled to 1), so the molecules move almost inde- 
pendently. Hence neglecting the elastic term in the Hamil- 
tonian, the partition sum in this regime is asymptotically 
given by 



l=1 Jo V 
which leads to the mean separation 



1 



(4.20) 



(4.21) 



In the same way, one obtains the surface roughness £j_ ~ 
1/Q. These results are displayed in Fig. |4.4| . 

Now, let u s in spect the probability distribution V(l), 
as defined in (4.9), in the SOS model. Except for its tails 
far away from the hard wall, V(l) is rather well fitted by 
a Gaufiian distribution of the form 



x 2< 




■ exp 



2£i 



for I > 
for I < 0. 



(4.22) 



This is shown in Fig. 4.5, where the rescaled func 



tions £,±V((l — l m ax)/£,±) collapse almost onto a single 
curve which coincides well with the Gaufiian distribution 
Vfi%(l). For distances I < £, i.e., on the side facing the hard 
wall, V(l) does indeed not deviate substantially from the 
Gaufiian. Hence, our argument leading to the prediction 
1/2 = 1 for the contact exponent as in (4.12) also applies to 
the SOS model. The scaling behavior for the probability 
for local pair contacts "P26 is fully confirmed by numerical 
results which have been obtained over a pressure range of 
five orders of magnitude, see Fig. L6. 
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Fig. 4.5. Data collapse £±V((l — lmax)/£±) for two surfaces in 
the SOS model, for 11 different values of the external pressure 
in the range Q — 5 x 1CP 5 — 7 x 10~ 2 . Imax is the most probable 
distance between the two surfaces. The topmost curve, which 
most strongly deviates from the Gaufiian, corresponds to the 
highest pressure value. 
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Fig. 4.6. Contact probability Vib versus pressure Q for the 
SOS-model. The dotted line corresponds to the value V2 = 1 
for the contact exponent. It turns out that Vib = cxQ with 
a ~ 1. 



For I > £, on the other hand, substantial deviations 
from the Gaufiian behavior are observed since the his- 



1 surface + hard wall 


.<&'' 


2 surfaces 




P 2b = Q 







»'' 
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tograms have an exponential tail; see also Fig. 4.7. This 
tail results from isolated protrusions, which can be under- 
stood as follows. The energy of an isolated protrusion of 
height I is given by E(l) = J2n.n. \h ~ l\/(n n Jsc), which, 
if all neighbouring sites i have heights ij < I, simplifies 
to E(l) = \l nn - l\/l sc , with l nn = n~l J2 n .n. l i- Together 
with the external pressure term Q, the probability of a 
single prot rusi on is hence proportional to e~ z ( 1+ ®\ The 
fits in Fig. 4.7 confirm this simple explanation. 
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Fig. 4.7. Probability distribution V{z) for one fluctuating 
surface with hard wall interaction within the SOS model. The 
tails have been fitted by exponentials of the form ce _z( - 1+<5 ', 
where the amplitude c represents a fit parameter. 



length lh y , which introduces a new length scale into the 
problem. 

Linear functional renormalization of the hydration in- 
teraction shows that the interplay between the protrusion 
length l pr and the hydration length lh v leads to two dif- 
ferent regimes Q: For l pr > 2l} ly and l pr < 2lh y , one has 
a protrusion regime and a hydration regime, respectively. 
In the protrusion regime, the asymptotic behavior is the 
same as in the case of a hard-wall interaction alone. In 
the hydration regime, on the other hand, one finds quali- 
tatively new behavior. 

Within the hydration regime, the relation between the 
pressure P and the interfacial roughness £j_ has again the 
Gaufiian form as given by ( 4^5 ) but with the new pressure 
amplitude P2± = lh y S/a 2 . The relation between the pres- 
sure and the mean distance £, on the other hand, is now 
given by 

P w P 2 e~ e/h (5.2) 



which is governed by the new length scale 

k = ky [1 + (Ipr/Zhy) 2 ] ■ (5.3) 

The pressure amplitude Pi has the parameter dependence 
P 2 = (Z/ay Vfl~ p l 2 ^ 1 (5.4) 



5 Two surfaces with soft wall interactions 

Now, we will study two surfaces which are governed by 
the potential 

V{1) = V hw {l) + V hy exp[-l/l hy ] + PI (5.1) 

which includes the hydration interaction Vh y e~ l / lhy . The 
decay length of this soft wall is given by the hydration 



with p = 1 - Ihy/h- 

Note that the length scale l± as given by (|5.3| ) depends 
both on the protrusion length l pr and on the hydration 
length lh y . For lh y ^> l pr , one has It ~ l hy ; i.e., in the 
limit in which the hydration potential has a much longer 
range than the effective protrusion potential, the mean 
separation is solely determined by the balance between 
the applied pressure and the repulsive hydration forces. 
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Fig. 5.1. Discrete Gaufiian model: Mean distance (z) versus 
pressure Q for several hydration lengths Zh y for one surface 
fluctuating against a soft wall, with a linear dimension Ln = 64 
and potential strength Uhy ~ 5.7. The data are fitted by the 
functional form z = z t ln(Q2/Q). 




Fig. 5.2. Discrete Gaufiian model: Rough ness £^ versus pres- 
sure Q from the same MC data as in Fig. 5.1 For the fit, the 
value Zp r /2 = 0.15 has been used, corresponding to z pr ~ 0.55. 
Note the finite size effects for Q < 0.02. 



5.1 GauBian models 

Let us introduce two additional dimensionless quantities 
for the soft wall interaction via 



Zhy — lhy /Is 



and 



U hy = AoVhy/T. (5.5) 



The above results for the pressure-dependence of t and £j_ 
are confirmed by MC si mula tions of the discrete Gaufiian 



model, as shown in Fig. 5.1 and Fig. 5.2 



As we have already shown in the preceding section, 
the probability distribution V(l) for the discrete Gaufiian 
model with a hard wall interaction has no exponential tail. 
In the presence of the hydration interaction, this distribu- 
tion is well fitted by the functional form 



Vcauil) ~ exp 



2£i 



U hy e 



(5.6) 



where the influence of the direct hydration potential has 
been taken into account in the simplest possible manner. 
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Fig. 5.3. Typical histograms for the Gaufiian model with 
a soft wall. Sample size 64x64, Uhy = 5.7, zh v = 1.41. The 
fitting curves have the functional form exp( — (z — (z)) 2 /2£± — 
Uh y e~ z / Zhy ). Note the absence of exponential tails. 



Note that the term that explicitely depends on the hy- 
dration potential in ( |5.6D is appreciable only for I < lh y , 
i.e. close to t he h ard wall, while I diverges in the limit of 
small P. Fig. 5.3 shows some fits of this form. Obviously, 
the tails have pronounced Gaufiian characterists, rather 
than exponential tails, on both sides. On the side facing the 
wall, the fits are excellent, while they are slightly poorer 
on the other side. 

An evaluation of the Gaufiian distribution Vcauil) at 
the position I = of the hard wall, leads to the contact 
probability 



-p 2b ^ e - f2 ml _ p(h/i pr ) 2 



(5.7) 



in the limit of small pressure P, w here the p ress ure de- 
pendence of £j_ and I as given by (4.8) and fl5.2j), resp., 
has been used, ignoring confluent logarithmic terms. 
Thus, the contact exponent i>i is found to be 



V-2 = {h/lprf 



(5.8) 



where the length scale It is given by ( |5.3| ). Note that for 
lhy = lpr/2, which defines the border line between the 
protrusion and the hydration regime, one has l t = l pr and 
hence f2 = 1 as in the protrusion regime. Therefore, the 
contact exponent v% is a continuous function of lhy 

The corresponding data from MC si mula tions of the 

An inspec- 



5.4 



discrete Gaufiian model are shown in Fig. 
tion of this figure reveals consi dera ble deviations from the 
predicted values as given by (5T8|). These deviations can 



be explained in the following way. 

First of all, the accessible regime of pressure values Q 
is quite restricted, covering only about one order of mag- 
nitude: For large Q > 1, we are deep in the mean field 
regime where the collective behavior is dominated by the 
external pressure. In the latter case, the molecules move 
almost independently and the probability distibution is 
~ Qexp(— Qz). This implies that the contact probability 
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Fig. 5.4. Contact exponents v-z as a function of the hydration 
length Zhy for the discrete Gaufiian model. In the inset, fits to 
MC data for the contact probability Vib versus pressure Q are 
shown. 
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Q 



Fig. 5.5. Contact exponents for the discrete Gaufiian model 
with soft wall interaction. The upper curve corresponds to 
Zhy = 1.4, the lower one t o Zh v = 2.1. For a detailed discussion 
of this figure, see section 5.1 



scales linearly with Q, not depending on the details of the 
hydration force. For small Q < 0.01, on the other hand, 
the finite siz e of the sample effects the interface roughness 
see Fig. 5^2. This has a large effect on the histogram, 
making it 'narrower' than it would be for larger systems, 
while the mean distance £ is essentially unaffected by the 
finite size. This lowers the observed contact probability 
Vib ~ ^(0). For this reason, we expect to extract a con- 
tact exponent from the data which is larger than the true 
asymptotic value. This behavior is expected for all values 
of Zhy, hence all data will tend to lie above the theoretical 
line. 

Secondly, the asymptotic form for V 2 b as given by fl5.7| ) 
is difficult to extract from intermediate values of Q be- 
cause of large corrections to the leading scaling behav- 
ior which arise as long as the pressure amplitudes Q 2 
and Qil take different values. Indeed, one has V 2 b ~ 
exp[— £ 2 /2£,p erp ) with the asymptotic behavior 



2^3 



;2 
_L 

P 

pr 

H 
I 2 

pr 



In 2 (Q 2 /Q) 
HQ 2 ±/Q) 



(5.9) 



In 



Ql 
Q2±Q 



In 2 (Q 2 /Q 



2-U 



InQ 



O 



ln 2 Q 



The leading order correction vanishes only for Q 2 = Q 2 ± ■ 
Using the values for Q 2 and Q 2 ± as measured in the MC 
simulations of the discrete Gaufiian model, we find that 
these corrections affect our data both for Zh y = 1.4 and 
for Zhy = 2.1. 

The measured contact exponents are however consis- 
tent with the actual values of exp[— £ 2 /2^ 2 _^] as determined 
from the measured values of i and In order to show 
this, we have plotted exp(-Zj ln 2 (Q 2 /Q) / I 2 ln(Q 2 ±/Q)), 



corresponding to the first equality in fl5.9| ), as a function 
of Q, see Fig. |5j| (full lines). The values for Q 2 and Q 2 ± 
were taken from the MC data. For comparison, we have 



included the expected asymptotic scaling function V 2 b ~ 
Q z t /V, with the contact exponents v 2 = z 2 / z 2 r = 7.14 
for Zhy — 1.4 and v 2 —\hA for Zh y =2.1 (dashed lines). Fi- 
nally, for the latter value Zh y =2.1, the apparent scaling 
in the intermediate regime 0.03 < Q < 0.2 with the expo- 
nent v 2 = 13 (as determined from MC data, see Fig. 5.4) 
is displayed (dashed- dotted line). Obviously, in the pres- 
sure regime that was accessible to MC simulations, this 
effective scaling exponent coincides well with the theoret- 
ical prediction. Note that this effect leads to data that are 
lower than the theoretical value for large values of Zhy 

Together, these two effects explain the deviations of 
the measured values for v 2 from the theoretical prediction 
v 2 = (k/lpr) 2 - 



5.2 Solid-on-solid models 



The SOS model with soft wall interactions has been ex- 



pressed in terms of the dimensionless variables (5.5) and 
(4.16). Using this dimensionless form of the model, we 
have first studied the critical behavior of the mean separa- 
tion I and the roughness £j_ . The corres pon ding MC data 
for Zhy > Zp r /2 are displayed in Fig. |5.6| and Fig. |5.7| . 
In these figures, we compare the critical behavior of (I) 
one fluctuating surface interacting with a hard wall and 
of (ii) two fluctuating surfaces. Inspection of these figures 
shows that both cases exhibit analogous critical behavior. 
We have also studied the critical behavior of I and £j_ for 
smaller values of zy ty . All of these data are again consistent 
with two distinct scaling regimes: a hydration regime for 
Zhy > Zp r /2 and a protrusion regime for Zh y < z pr /2. 

Next, let us consider the contact probability V 2 b- In 
the protrusion regime, the mean separation between the 
membranes and their width are determined solely by the 
protrusion length l pri so the hydration interaction should 
not alter the pressure dependence of 7^26 • Thus, one ex- 
pects again v 2 = 1 . This is confirmed by numerical results 
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5.6. Hydration regime in the SOS model: The length 
z t is plotted versus z pr /zh y , for one and for two flue- 



Fig. 

scale Zt is plotted versus z pr /Zh y , tor one 
tuating surfaces. The broken curve corresponds 
diction of functional renormalization as given by (5.3) 
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Fig. 5.8. Protrusion regime in the SOS model: Contact prob- 
ability Vib versus pressure Q with Zh y < z pr /2. The data con- 
+ Zpr/izhy), with z pr = 1.54 for one surface and nrm ^ ne universal value U2 = 1 for the contact exponent in this 



= 1.07 for two surfaces. 
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Fig. 5.7. The pressure amplitudes Q2 and Q2X versus hydra- 
tion length z^, for the SOS model of one and two fluctuat- 
ing surfaces, respectively. For comparison, the formulae from 
Eq. are plotted. 



that give the values — 1.015 and 1.025 for zty = 0.3 
and 0.5, respectively, which coincides with t he e xpected 
value within the numerical precision, see Fig. |5.8| . 

The situation changes when one enters the hydration 
regime with z^y > z pr /2. Now, the mean distance is pro- 
portional to It, while the roughness is still determined by 
l pr . In the preceding chapter, we have found the value 
v 2 — (h/lpr) 2 for the hydration regime of the Gaufiian 
model. However, we will now argue that, in the hydration 
regime of the SOS model, the contact probability V-x, is 
strongly affected by exponential tails of the probability 
distribution arising from molecular protrusions. 

In the protrusion regime of the SOS model, the proba- 
bility distribution V(l) has a pronounced exponential tail 
on the side of the interface which does not face the wall, 
while on the side facing the wall, it is essentially Gaufiian 
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Fig. 5.9. Typical probability distributions V(z) for the SOS 
model with one fluctuating surface and a soft wall. The sam- 
ple size was 64x64 sites, and the hydration interaction pa- 
rameters were set to Uhy = 2.0 and Zhy = 4.0. The tails are 
well represented by a mean field like probability distribution 
V(z) ~ exp(-|z— (z} \ — Qz- U hy e~ z/zh -y). Each tail has been 
fitted individually. In each fit, the only fit parameter was a con- 
stant factor. The difference to the histograms for the Gaufiian 
model, where the exponential tails are absent, is obvious. 



shaped. This is the same behavior as observed in the case 



of pure hard wall interactions, as discussed in section 4.2 
Now, in the hydration regime of the SOS model, expo- 
nential tails a re fo und on both sides of the histogram, as 
shown in Fig. 5.9. This observation inspires a theoretical 
prediction for the contact exponent that yields a much 
better description of the data, based on the assumption 
that these tails are the result of isolated protrusions which 
can be considered as 'rare' events. 

Thus, we should distinguish two different types of ex- 
citations: (i), Collective fluctuation modes that lead to a 
Gaufiian shaped probability distribut ion and, consequently, 
to the contact exponent as given by (5.8); and (ii) Isolated 
protrusions of single molecules which are absent in the dis- 
crete Gaufiian model. 
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Now, let us assume that the tails of the distribution 
V(l) and hence the contact probability Vib is dominated 
by these rare protrusions. The simplest estimate for Vib 
is now given by 



7 SOS 
2b 



exp[-£/l s 



(5.10) 



which describes the probability for a single protrusion 
mode to extend over the distance I between the inter face 
and the hard wall. With the scaling of t as given by ( |5.2| ) 
and (5.3), we now find the contact exponent 



Vl — h/lsc — {lhy/lsc)[l + (Ipr/^lhy)" 2 



(5.11) 



and 



Note that this value depends on the two ratios lhy /Is 

Ihy/lpr • 

Thus, the GauBian fluctuations and the single pro- 
trusion modes lead to a different asymptotic behavior of 
the contact probability Vib- The true asymptotic behav- 
ior will be determined by those fluctuations which dom- 
inate Vib- This competition leads to another characteris- 
tic length scale ij£ which separates two different scaling 
regimes within the hydration regime. 

Direct com paris on o f the contact probabilities Vib as 
given by (5.7) and ( 5.10 ) leads to the characteristic length 
scale 

/ "' ' ~ ^ (5.12) 



/* — 

L hy ~ 



T 



'pr 



i 



With z pr = 1.54, this implies l* hy ~ 2.71 x l pr /2 ~ 2.09 l sc . 

The GauBian fluctuations and the single protrusion 
modes dominate for lh y < lt y and h iy > It , respectively. 



l hy <■ L hy dna ''hy 

This leads to the contact exponent 



V2 = {k/lpr) 2 for l pr /2 < l h y < l* hy 
= k/lsc for l h y > lhy, 



In Fig. 5.10 , MC data for the contact exponent in the 
SOS mod el ar c shown, together with the theoretical pre- 
dictions (5.13). The different values of the contact expo- 



nent v<x are summarized in Table 1. 
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Fig. 5.10. MC data for a system with one flexible surface 
and a soft wall. The numerical data for the contact exponent i>2 
are compared with i/ 2 = 1 in th e protrusion regime, and with 
the functions as given by (5.13) in the hydration regime: The 
dotted line corresponds to the prediction from GauBian fluc- 
tuations, U2 — zt/zp r , and the dashed-dotted line to V2 = Zt. 
Here, we have used z pr = 1.54. Note that these two curves cross 
at z^y = 2.09. Inset: The MC data for P26 are plotted versus 
the scaled pressure Q and fitted against Vib ~ Q" 2 • Apart from 
Zhy = 8.0, all data cover more than an order of magnitude in 
the pressure range and several orders of magnitude in Tib- 



Table 1. The contact exponent ^2 for the discrete GauBian 
and the solid-on-solid models. 
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